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Claims 

[d] What is claimed is: 

l.A method for characterizing a subsurface formation, 
comprising the steps of: 

disposing within a borehole a logging instrument 
equipped with at least first transmitter and receiver an- 
tennas spaced apart by a first distance, at least one of 
the first antennas having a tilted magnetic dipole with 
respect to the longitudinal axis of the instrument, the 
antennas being oriented about the axis of the logging 
instrument such that the at least one tilted magnetic 
dipole corresponds to a first azimuthal angle; 
azimuthally-rotating the logging instrument within the 
borehole; 

while the logging instrument is rotating, activating the 
first transmitter antenna to transmit electromagnetic en- 
ergy into the formation; 

while the logging instrument is rotating, directionally 
measuring the first voltage signals associated with the 
transmitted electromagnetic energy using the first re- 
ceiver antenna, as a function of the azimuthal orientation 
of the logging instrument, so as to determine the az- 
imuthal variation of the measured first voltage signals; 



and 

fitting the azimuthal variation of the measured first volt- 
age signals to approximate functions. 

[c2] 2.The method of claim 1, wherein the fitting step is exe- 
cuted while the first voltage signals are being measured. 

[c3] 3.The method of claim 2, further comprising the step of 
stopping the fitting when convergence has been 
achieved. 

[c4] 4.The method of claim 3, wherein the activating, mea- 
suring, and fitting steps are repeated to execute subse- 
quent acquisition cycles. 

[c5] 5.The method of claim 1, wherein the fitting functions 
are sinusoids defined by coupling components of the 
first transmitter antenna"s magnetic dipole and first re- 
ceiver antenna"s orientation vectors. 

[c6] 6.The method of claim 5, wherein the coefficients of the 
fitting components are functions of earth formation pa- 
rameters including at least one of resistivity of formation 
beds, location of the logging instrument, borehole devi- 
ation, azimuth angle at the location of the logging in- 
strument, and a combination thereof. 

[c7] 7.The method of claim 6, wherein the fitting coefficients 



include constant, sinO, cosO, sin20 and cos20 terms that 
define an iterative fitting algorithm useful for determin- 
ing the azimuthal dependence of the directional mea- 
surements. 

[c8] 8.The method of claim 7, wherein the iterative fitting al- 
gorithm is used for selected real-time directional mea- 
surements having utility in geosteering. 

[c9] 9.The method of claim 1, where the fitting coefficients 
are determined using a Fast Fourier Transform. 

[dO] lfj.The method of claim 1, wherein the logging instru- 
ment is further equipped with second transmitter an- 
tenna and receiver antennas spaced apart by the first 
distance, the second transmitter having a magnetic 
dipole whose tilt corresponds to the tilt of the first re- 
ceiver antenna and the second receiver antenna having a 
magnetic dipole whose tilt corresponds to the tilt of the 
first transmitter antenna such that at least one of the 
second antennas has a tilted magnetic dipole, the second 
transmitter and receiver antennas being oriented about 
the axis of the logging instrument such that the at least 
one tilted magnetic dipole corresponds to a second az- 
imuthal angle, and further comprising the steps of: 
while the logging instrument is rotating, activating the 
second transmitter antenna to transmit electromagnetic 



energy into the formation; 

while the logging instrument is rotating, directionally 
measuring the second voltage signals associated with 
the transmitted electromagnetic energy using the second 
receiver antenna, as a function of the azimuthal orienta- 
tion of the logging instrument, so as to determine the 
azimuthal variation of the measured second voltage sig- 
nals; 

fitting the azimuthal variation of the measured second 
voltage signals from the second receiver to approximate 
functions. 

[c11] ll.The method of claim 10, wherein the second az- 
imuthal angle differs from the first azimuthal angle by 
substantially 90 degrees. 

[d2] i2.The method of claim 10, wherein the second az- 
imuthal angle is substantially equal to the first azimuthal 
angle. 

[d3] l3.The method of claim 10, wherein the fitting step is 
executed while the second voltage signals are being 
measured. 

[d4] i4.The method of claim 13, further comprising the step 
of stopping the fitting when convergence criterion has 
been achieved. 



[d5] l5.The method of claim 14, wherein the activating, mea- 
suring, and fitting steps are repeated to execute subse- 
quent acquisition cycles. 

[d6] l6.The method of claim 10, wherein the fitting functions 
are sinusoids defined by coupling components of the 
first transmitter antenna"s magnetic dipoles and first re- 
ceiver antenna"s orientation vectors, and by coupling 
components of the second transmitter antenna"s mag- 
netic dipoles and second receiver antenna"s orientation 
vectors. 

[d7] l7.The method of claim 16, wherein the coefficients of 
the fitting components are functions of earth formation 
parameters such as resistivity of formation beds, loca- 
tion of the logging instrument, borehole deviation, and 
azimuth angle at the location of the logging instrument. 

[d8] l8.The method of claim 17, wherein the fitting coeffi- 
cients include constant, sinO, cosO, sin20 and cos20 
terms that define an iterative fitting algorithm useful for 
determining the azimuthal dependence of the directional 
measurements. 

[d9] lg.The method of claim 18, wherein: 

the measured first and second voltage signals are com- 
plex voltage signals; and 



further comprising the steps of: 

calculating the phase-shift and attenuation values from 
the fitting coefficients for the measured first and second 
voltage signals; 

combining the calculated phase-shift and attenuation 
values for the measured first and second voltage signals 
to generate a symmetrized or anti-symmetrized mea- 
surement. 

[c20] 20.The method of claim 19, wherein the combination of 
calculated phase-shift and attenuation values is achieved 
by one of a sum and a difference. 

[c21] 21.The method of claim 18, wherein: 

only one of the first antennas has a tilted magnetic 
dipole; and 

further comprising the step of: 

characterizing the noise of the measured first and sec- 
ond voltage signals using the second harmonic coeffi- 
cients. 

[c22] 22.The method of claim 18, wherein: 

each of the first antennas has one of a tilted and a trans- 
verse magnetic dipole; and 

the fitting coefficients include third harmonic coeffi- 
cients; and 

further comprising the step of: 



characterizing the noise of the measured first and sec- 
ond voltage signals using the third harmonic coeffi- 
cients. 

[c23] 23.The method of claim 18, further comprising the step 
of characterizing the noise of the first and second mea- 
sured voltage signals by combining the first and second 
measured voltage signals. 

[c24] 24.The method of claim 19, wherein the phase-shift and 
attenuation values are obtained by taking the logarithm 
of the ratio of the complex voltage signals obtained from 
the fitting expression at two azimuthal angles. 

[c25] 25. The method of claim 24, wherein the two azimuthal 
angles are 0 and 180 degrees from a determined bed- 
ding azimuth. 

[c26] 26.The method of claim 18, wherein the second anten- 
nas are symmetric with respect to the first antennas, and 
further comprising the steps of: 
suspending rotation of the logging instrument; 
determining the azimuth of a bed of interest by combin- 
ing the first and second antenna couplings; 
determining the constant and first harmonic coefficients 
from measured first and second voltage signals acquired 
when the instrument is not rotating; 



using the constant coefficients to execute the fitting step 
when the logging instrument is rotating. 

[c27] 27.The method of claim 26, further comprising the steps 
of: 

updating the azimuthal variations of the measured first 
and second voltage signals; and 

re-calculating equivalent voltages when the tool is in the 
plane of the bedding. 

[c28] 28.The method of claim 19, wherein the phase-shift and 
attenuation values are obtained from log ((cO + cl)/(cO - 
cl)), where cO and cl are the complex fitting coefficients 
taking a determined bedding azimuth as an angle refer- 
ence. 

[c29] 29.The method of claim 7, wherein the iterative fitting 
algorithm includes the steps of: 



initialize F„ and U^, 

for m = 1 to Nsamplss 

next m; 
retum(Uy, 

Ifyamples is the total number o f samples acquired in one cycle, 

is the dimension of the approximate function vector (number of approximation 
functions), ^ 
U is the vector of fitting coefficients of dimension , 

7 is the vector of approximate function values at each measure position of dimension 
and M ■■■■■ M 

P is a matrix of dimension 



[c30] 30.The method of claim 29, wherein the iterative fitting 
algorithm determines if the fit error is below a pre- 
defined threshold, and if U converges to a value that is 
representative of the fitting coefficients. 

[c31] 31.The method of claim 30, wherein the iterative fitting 
algorithm employs an integer implementation. 

[c32] 32.The method of claim 29, further comprising the step 
of using the fitting coefficients to determine the orienta- 
tion of a formation bed. 

[c33] 33.The method of claim 32, wherein: 

the measured first and second voltage signals are com- 
plex voltage signals; and 



the orientation of the formation bed with respect to the 
azimuthal angle reference for each channel of directional 
measurement is determined according to: 



= tan _1 [ 



where Cjjisthe real or imaginary part of the coefficient of 5in$, and C; f is the coefficient 
of cos ^ from the fitting. 



[c34] 34. The method of claim 33, further comprising the step 
of calculating a common azimuthal angle for the first 
and second voltage signals using weighted averaging of 
the fitting coefficients for real and imaginary parts of the 
measured voltage signals. 

[c35] 35.The method of claim 34, further comprising the step 
of calculating the amplitude and phase of the measured 
voltage signal at an assumed normal direction to a bed 
boundary of interest. 

[c36] 36.The method of claim 35, further comprising the step 
of determining phase shift and attenuation by taking 
propagation measurements for two azimuth angles. 

[c37] 37.The method of claim 36, further comprising the ap- 
plication of an inversion technique to interpret the direc- 
tional measurements. 



[c38] 38.The method of claim 36, wherein the azimuth angles 
are 

[c39] 39.The method of claim 36, further comprising the step 
of combining the signals from the fitting coefficients for 
the first and second voltage measurements, to produce 
signals necessary to determine the distance to bed 
boundaries of interest. 

[c40] 40.A method for characterizing a subsurface formation, 
comprising the step of: 

cross-plotting two directional logging measurements ac- 
quired from an instrument disposed in a borehole inter- 
secting the formation to obtain a distance to at least one 
formation boundary and a resistivity for at least one for- 
mation bed. 

[c41] 4i. The method of claim 40, wherein: 

the cross-plotting is achieved using a one-boundary 
model; 

the obtained resistivity is the shoulder-bed resistivity; 
and 

the obtained distance is the closest distance to the 
shoulder-bed. 

[c42] 42.The method of claim 40, wherein the cross-plotting 



step includes the steps of: 

defining an appropriate model; 

selecting appropriate directional measurements; 

inputting the selected measurements to the defined 

model to generate the cross-plot; and 

generating a visual representation of the cross-plot. 

[c43] 43.The method of claim 42, further comprising the step 
of updating the cross-plot with further real-time mea- 
surements. 

[c44] 44.The method of claim 43, further comprising the step 
of: 

using cross-plot definitions and the real-time measure- 
ments to determine a resistivity for at least one forma- 
tion bed and a distance to at least one formation bound- 
ary. 

[c45] 45.The method of claim 40, wherein the obtained dis- 
tance and resistivity are used to make drilling decisions. 

[c46] 46.A method for characterizing a subsurface formation, 
comprising the step of: 

cross-plotting a resistivity and a directional measure- 
ment determined using an instrument disposed in a 
borehole intersecting the formation to obtain a distance 
to at least one formation boundary and a resistivity for at 



least one formation bed. 

[c47] 47. The method of claim 46, wherein: 

the cross-plotting is achieved using a one-boundary 
model; 

the obtained resistivity is the shoulder-bed resistivity; 
and 

the obtained distance is the closest distance to the 
shoulder-bed. 

[c48] 48.A method for characterizing a subsurface formation, 
comprising the step of: 

cross-plotting a resistivity and two directional measure- 
ments determined using an instrument disposed in a 
borehole intersecting the formation to obtain a distance 
to at least one formation boundary and a resistivity for at 
least two formation beds. 

[c49] 49. The method of claim 48, wherein: 

the cross-plotting is achieved using a one-boundary 
model; 

the obtained resistivities are the bed and shoulder-bed 
resistivities; and 

the obtained distance is the closest distance to the 
shoulder-bed. 

[c50] 50.The method of claim 48, further comprising the steps 



of: 

selecting an appropriate inversion model for the selected 
real-time directional measurements; 
verifying that the selected model is consistent with other 
information; and 

using the verified model to make drilling decisions. 

[c51] 51.The method of claim 50, wherein the model-selection 
step includes running at least one model-type selected 
from the set of: 

homogenous isotropic (single parameter: resistivity); 
homogenous anisotropic (two parameters: Rh and Rv); 
single boundary isotropic formation, boundary above or 
below (three parameters: Rbed, Rshoulder and distance 
to boundary); 

single boundary anisotropic formation, boundary above 
or below (four parameters: Rbed_h, Rbed_v, Rshoulder 
and distance to boundary); 

two boundary isotropic formation three parameters: (five 
parameters: Rbed, Rshoulder_up, Rshoulder.down and 
distance to boundary above and below the tool); and 
two boundary anisotropic formation three parameters: 
(six parameters: Rbed_h, Rbed_v, Rshoulder_up, Rshoul- 
der.down and distance to boundary above and below the 
tool). 



[c52] 52.The method of claim 50, wherein the model selection 
step includes creating a visualization of the selected di- 
rectional measurements. 

[c53] 53.The method of claim 50, wherein the model selection 
step includes use Akaike Information Criterion to penal- 
ize the model complexity. 

[c54] 54.The method of claim 50, wherein the model selection 
step includes: 

identifying known formation parameters; 
interactively choosing the models with which to invert 
the selected directional measurements; and 
selecting the simplest model that fits the known infor- 
mation. 

[c55] 55.The method of claim 50, wherein the verifying step 
includes the step of: 

comparing the selected model to known geological char- 
acteristics and other measured formation parameters. 

[c56] 56.The method of claim 55, wherein the verifying step 
further includes the step of updating the selected model 
if the selected model is not consistent with the known 
information. 

[c57] 57.The method of claim 56, wherein the updating step 
includes the steps of: 



refining the selected model based upon one of trends, 
prior knowledge, external information, and a combina- 
tion thereof; 

selecting parameters to be inverted; 

defining ranges for all parameters to be inverted; 

updating the inversion model by adding more formation 

beds; 

re-weighting or eliminating some of the selected real- 
time directional measurements; and 
re-inverting the resulting real-time directional measure- 
ments to the updated model. 

[c58] 58.An apparatus for measuring characteristics of earth 
formations surrounding a borehole, comprising: 
a logging instrument adapted for disposal within the 
borehole, the logging instrument having a longitudinal 
axis and being equipped with first and second transmit- 
ter-receiver antenna pairs; 

the first transmitter-receiver antenna pair comprising 
a first transmitter antenna having a magnetic dipole ori- 
ented in a first direction with respect to the longitudinal 
axis of the logging instrument, 

a first receiver antenna located a first distance away from 
the first transmitter antenna and having a magnetic 
dipole oriented in a second direction, the first and sec- 
ond directions being different, 



the magnetic dipoles of the first transmitter and receiver 
antennas defining a plane that includes the longitudinal 
axis of the logging instrument, 

the second transmitter-receiver antenna pair comprising 
a second transmitter antenna having a magnetic dipole 
oriented in the second direction with respect to the lon- 
gitudinal axis of the logging instrument, 
a second receiver antenna located the first distance away 
from the second transmitter antenna and having a mag- 
netic dipole oriented in the first direction, 
the magnetic dipoles of the second transmitter and re- 
ceiver antennas defining a plane that includes the longi- 
tudinal axis of the logging instrument, 
a toolface sensor for continuously indicating the az- 
imuthal orientation of the logging instrument; and 
a controller for controlling the first and second transmit- 
ter-receiver antennas pairs so as to selectively transmit 
electromagnetic energy into the formation and measure 
the voltage signals associated with the transmitted elec- 
tromagnetic energy as a function of the azimuthal orien- 
tation of the logging instrument. 

[c59] 59.The apparatus of claim 58, wherein the second trans- 
mitter-receiver antenna pairs are oriented at a first az- 
imuthal angle with respect to the first transmitter-re- 
ceiver antenna pairs about the longitudinal axis of the 



logging instrument. 

[c60] 60.The apparatus of claim 58, further comprising: 

a CPU for processing the measured voltage signals within 
the borehole; 

a telemetry apparatus for transmitting the measured sig- 
nals and CPU-processed results from the borehole to the 
surface; and 

a surface system for further processing measured signals 
together with other measurements to generate and dis- 
play selected parameters of a consistent earth model. 

[c61] 61.The apparatus of claim 58, wherein the formation 
characteristic is resistivity. 

[c62] 62.The apparatus of claim 58, wherein the formation 
characteristic is the geometry information of the earth, 
including the dip, azimuth, and bed thickness. 

[c63] 63.The apparatus of claim 58, wherein the first direction 
is substantially collinear with the longitudinal axis of the 
logging instrument. 

[c64] 64.The apparatus of claim 58, wherein the second direc- 
tion is substantially collinear with the longitudinal axis of 
the logging instrument. 

[c65] 65.The apparatus of claim 58, wherein the first direction 



is substantially 45 degrees from the longitudinal axis of 
the logging instrument. 



[c66] 66.The apparatus of claim 58, wherein the second direc- 
tion is substantially 45 degrees from the longitudinal 
axis of the logging instrument. 

[c67] 67.The apparatus of claim 58, wherein each of the trans- 
mitters and receivers have transceiver capabilities. 

[c68] 68.The apparatus of claim 59, wherein the first az- 
imuthal angle is substantially 90 degrees. 

[c69] 69.The apparatus of claim 58, wherein the first and sec- 
ond transmitter-receiver antenna pairs are located at the 
same physical positions on the logging instrument. 

[c70] 70.The apparatus of claim 58, wherein the toolface sen- 
sor employs magnetometers to indicate the azimuthal 
orientation of the logging instrument with respect to 
earth"s magnetic north. 

[c71] 71.The apparatus of claim 58, where the toolface sensor 
employs gravitation sensors to indicate the azimuthal 
orientation of the logging instrument with respect to the 
earth"s gravity vector. 



